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A one dimensional trap model for a thermally activated classical particle is introduced to si- 
mulate driven dynamics in presence of "ageing" effects. The depth of each trap increases with the 
time elapsed since the particle has fallen into it. The consequences of this dynamical pinning are 
studied, and velocity-force characteristics are numerically obtained. A special attention is paid to 
the situation where the particle is pulled with a spring to ensure a finite average velocity. In the 
low velocity regime, the presence of a broad distribution of trapping times leads to suppression of 
hnear response, replaced by a threshold or by sublinear dynamics. A regime of strong fluctuations 
is obtained when the particle is driven at intermediate velocities. 



I. INTRODUCTION 

The out-of-equilibrium dynamics of glassy systems is 
a widely open subject. An important topic concerns 
models for which a spontaneous ageing behavior com- 
petes with an external field. This article aims to study 
a realization of the above paradigm in a simple "fric- 
tion" experiment where the two main features -ageing 
and external forcing- are present. For this purpose, a 
stochastic process describing a classical particle in a one- 
dimensional pinning potential is defined. The pinning 
potential consists of identical traps with time-dependent 
depths. The barriers between two neighbouring traps 
therefore increase with time in order to simulate age- 
ing effects, and the hopping rate between traps is based 
on Arrhenius dynamics. Some physical justifications for 
such time-dependent barriers are exposed in section ||. 
Indeed, by choosing adequately the time dependence of 
the trapping potential, one can account for anomalous 
slow diffusion situations related to broad algebraic dis- 
tributions of waiting times and leading to glassy behav- 
iors One can alternatively consider a simpler case 
where the barriers increases from a short-time value to a 
larger long-time value and discuss the related dynamical 
consequences. 

Our main purpose is to understand the driven dyna- 
mics of this system in the situation where the average 
velocity v is non zero. A quadratic potential moving at 
velocity v is therefore added, acting in a way similar to a 
spring pulling a body in a friction experiment. The mean 
spring extension defines the friction force F. The rela- 
tion between v and F is called the "Velocity-Force" char- 
acteristics, and is the central output for such a friction 
experiment. Other quantities of interest will be moni- 
tored like the distribution of waiting times between hops 
or the histogram of the "spring lengths" . In some situ- 
ations, the time-dependent pinning potential is found to 
enhance the fiuctuations of the particle position. 

The basic feature of this model consists in a competi- 
tion between the "ageing" pinning coming from the traps. 



and the renewal due to the driven dynamics. It results 
in a strongly non-linear v — F characteristics for a wide 
intermediate velocity regime. The more extreme case of 
logarithmic growing barriers leads to a non-ohmic be- 
haviour at low velocity, typical of glassy systems. In the 
latter case, a naive, "mean-field like" , regularization with 
a characteristic interruption time tc{v), fails in describ- 
ing the low-velocity properties of the model. The limit 
u — !• indeed appears to be singular and involves large 
pinning times ~ l/i;. 

Section II provides physical justifications of the model, 
detailed in Section |II]. A physical application is dis- 



cussed. Section [V focuses on the velocity-force charac- 



teristics induced by logarithmic growing barriers. Then, 
section studies some consequences of rapidly growing 
barriers on the particle's dynamics. 



II. A REVIEW OF MECHANISMS GENERATING 
AN EFFECTIVE TIME DEPENDENT PINNING 
POTENTIAL. 

A time dependent pinning potential. Many disordered 
systems like spin glasses, pinned random manifolds, dif- 
fusing particles in random media, involve a disordered 
energy landscape. In such systems, thermal effects, - 
diffusion and activation- compete against pinning effects. 
At low enough temperature, the out-of-equilibrium dy- 
namics is usually governed by activated barrier crossings 
in the configuration space. When the distribution of en- 
ergy barriers is broad, the time scale for reaching a ther- 
mal -Boltzmann- equilibrium may diverge with the size 
of the system. 

By following a particle during its out-of-equilibrium 
dynamics, one can attempt to define a function of time, 
with dimension of an energy, as the mean (or typical) 
height of the barrier that the system has to overcome, 
in order to expand further. This function "7i(i)" will be 
the starting point for this phcnomcnological approach of 
driven dynamics. 
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For instance, many aspects of the slow thermal diffu- 
sion of a particle in a one-dimensional quenched random 
force field, first introduced by Sinai' , have been partially 
reinterpreted as the consequence of such effective barri- 
ers, slowly increasing with time [Q. The random force 
generates a pinning potential -in 1 dimension- with un- 
bound local extrema, and the particle has to overcome 
higher and higher barriers in trying to reach its ther- 
mal equilibrium. Let i^, be the time interval elapsed 
since the beginning of the diffusion process, or waiting 
time. After t^, the particle is shown to be at equilibrium 
within a restricted area bounded by a barrier of height 
Ti.{tyj) ~ Tln{tyj) 1^,^. Such a growing effective barrier 
induces a strong dependence in i^,, or ageing, for the re- 
sponse and diffusion properties of the particle, together 
with an "anomalous" slow diffusion behavior. A gener- 
alization to other types of correlated disorder potentials 
leads to a variety of interesting regimes with sublinear 
characteristics 

In a related way, the ageing properties of a particle 
diffusing in a high dimensional phase space within traps 
with exponentially distributed depths, were extensively 
studied, in relation with the magnetic relaxation of spin 
glasses Again, after a waiting time t^,, the particle 
stays in a trap of typical depth Ti{tyj) ^ T\n{tw) 
Such a situation may occur for instance when describing 
the creep of an extended object like an elastic manifold. 
The pinning energy of this object increases as its internal 
degrees of freedom wander and find favorable configura- 
tions. One expects that the longer the object stays at the 
same place, the larger its pinning energy will be. A well 
known example involving pinning and motion of elastic 
lines concerns the fiux line problem in high Tc supercon- 
ductors ||]. 

However, we must keep in mind that many other sys- 
tems have slow relaxation properties which do not origi- 
nate from a barrier pinning mechanism. This is the case 
for a particle evolving in a high dimensional, random pin- 
ning potential |^ . In this system, thermal activated pro- 
cesses cannot explain the whole dynamics, because a non 
trivial relaxation subsists when the temperature goes to 
0. In such situations, the slow dynamics is attributed to 
"entropic effects" , i. e the system has more and more dif- 
ficulties to find favorable regions in its phase space, even 
if it is not separated from them by any energy barrier 
. The existence of growing effective energy barriers is 
therefore only a sufficient condition for ageing. 

In what follows, Section IV deals with particles dif- 
fusing among traps whose depth increases like Tln(i^), 
leading to a power law distribution for waiting times. In 
the last section , we focus on a less singular case where 
the time-dependent barrier Ti.(oo) < oo does not diverge. 

Detail of the friction experiment. In a friction exper- 
iment, one can either impose the external driving force 
F or the velocity drift v. In the first case, F is a con- 
stant and one measures the mean position {x{t)) of the 
particle. One can alternatively pull the particle with a 



spring, with the other end moving at constant velocity 
V. By doing so, one enforces only the average velocity, 
and the force is then related to the mean extension of the 
spring. 

Some systems make these two procedures be inequiva- 
lent. Some examples of random walks over random traps 
which exhibits only a sub-linear response for the mean 
displacement {x{t)) (xt"; a < 1 are reviewed in |T^. On 
the contrary, driving the particle with a spring ensures 
that the average displacement x{t) is linear in time, in 
other words the average velocity always exists, provided 
the average is taken on sufficiently long times. Never- 
theless, anomalous diffusion effects are still present and 
appears through the effective friction force, with the dis- 
appearance of the usual linear regime at low velocity. 
The study of such thermal "creep" for traps models gen- 
eralized in a way to allow the particle to be driven by a 
spring, will be the main purpose of section IV. As far as 
possible, results for the spring driven case will be com- 
pared to their "constant force" counterparts. 

On physical grounds, besides its "regularization" ef- 
fect, the quadratic confinement potential may refiect 
more complicated systems. For instance, when describ- 
ing a set of interacting particles, a crude mean-field 
elastic approximation produces a "cage potential" (self- 
averaging contribution of the whole system) which lim- 
its the excursion of a given particle This potential 
moves at velocity v with the center of gravity of the en- 
tire body. It can therefore be considered as a first step 
towards the incorporation of interactions in anomalous 
diffusion problems. 

As far as real "dry" friction experiments are concerned, 
it is known that phenomena similar to ageing occur at 
the contacts, and should be responsible for both the slow 
time dependence of threshold forces and for the "stick- 
slip" motion pj]Jl^ . Spring pulling and velocity depen- 
dent ageing are already present in Ref but their 
model is designed for a very different purpose. First, 
they introduce a very strong ageing leading to a reen- 
trant velocity-force characteristics and unstable motion. 
Then, they consider only a crude "mean-field" approach, 
called "adiabatic approximation" . 

In what follows we focus on overdamped, thermally ac- 
tivated motion, and especially on the low-velocity creep 
motion. In contrast with the model of Rcf jl^, the age- 
ing process is here described through a broad distribu- 
tion of trapping times, and the dynamics is simulated, 
giving access to nonlinear creep characteristics as well as 
to fluctuation effects. We provide details on the stochas- 
tic process and consider carefully the whole distribution 
of waiting times and spring lengths, which turn out to 
have a crucial importance for understanding the physical 
contents of this model: the dangers of a naive mean-field 
approach are emphasized in section IV. 
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III. TRAP DYNAMICS IN PRESENCE OF AN 
EXTERNAL DRIVING FORCE. 



LjQ exp 



exp(-/3H(0)) 



We consider a saw-tooth potential made of equidistant 
triangular wells (figure |l|) . The interval between two sites 
is constant, equal to d. The particle is supposed to hop 
towards its nearest neighbours sites, at a rate given by the 
Arrhenius law. Whereas the conventional "reaction rate 
theory" requires additional information about the spe- 
cific shape of the barrier, and provides corrective terms 
Ip^ , we will keep only the crude Arrhenius ratio. When 
a particle falls into a well, this one starts to increase its 
depth, in order to mimic an effective barrier grow- 
ing with time. The corresponding escape rate decreases 
accordingly, generating longer waiting times. A particle 
escaping from a well, immediately falls into one of its 
nearest neighbours. The particle is then always trapped, 
excepted during the negligible transit delay over the bar- 
rier. The saw-tooth potential is quite arbitrary, and the 
Arrhenius dynamics insensitive to this specific shape. 

An external potential -constant driving force or spring 
traction- tilts the saw-tooth potential, reduces one of the 
two barriers and increases the instantaneous escape rate. 
The barrier is strongly reduced by any external bias, un- 
til it vanishes at the critical force Fc{t) — 2TL{t)/d. Let 
h~{T) and h^ir) be the height of the left and right bar- 
rier (figure T is the time elapsed since the arrival 
of the particle in the occupied trap. By adding the left 
and right "channels" , the total escape rate wtot reads, {(3 
standing for the inverse temperature T~^, loq for a trial 
frequency) : 

wtotir) ^'^eM-Ph+{T)) + ^exp(-/?/i-(r)). (1) 

In absence of external force, the left and right barriers 
coincide and are set equal to T-L{t). By choosing 'H{t) ar- 
bitrarily, one can generate all the possible waiting time 
distributions. 71 constant corresponds to an ordinary ac- 
tivated diffusion over identical traps, leading to a diffu- 
sion constant D and mobility 

eM-Pn) , M = ^ eM-Pn). (2) 

An external force F changes H — > = H — Fd/2; 
h~ = Ti, + Fd/2}, leading to the constant force F escape 
rate: 



(3) 



A spring of stiffness fc, with a head moving at constant 
velocity v and the particle located at x{t) exerts a force 
F{t) = kl{t) with l{t) = vt — x{t). The total escape rate 
Wtotir, l(^r=o)) reads : 



wi exp[-/3W(T) + fSniO)] ■ ch 



f Pkdivt - x(t)) 
V 2 



It's essential to distinguish between the total time t 
and the delay r relative to the more recent jump of the 
particle, which occurred at t — r. Once the value of 1(t=q) 
is known, at the time of arrival of the particle, (^ deter- 
mines the probability distribution for the waiting time 
before the next jump. 

As expected, the external force affects the escape rate. 
At large times, the hyperbolic cosine dominates the ex- 
ponential term and make the particle eventually depin. 
Depinning is fast when Wtot becomes larger than 1, which 
happens when the force is greater than the critical ratio 
2niT)/d. 

Let's denote by n(t) the probability for the particle 
to stay within the well for a time r longer than t. n(i) 
decays exponentially as : 



Hit) = exp 



ds wtotis,l(s=o)) 



(5) 



By definition, 7r(i) — — dn(^)/dt is the probability dis- 
tribution of waiting times. The knowledge of n(t) allows 
to perform "molecular dynamics" simulations by follow- 
ing a procedure close to BKL In this simulation 
scheme, one computes at each step, the waiting time r 
before the next jump. This is achieved by converting 
a uniform random number r e [0,1] into T(r), i.e. by 
inverting n(T) — r. 

The stochastic process consists in a sequence of jumps, 
labeled by i. Each jump has a random direction 
Gi ~ + / — 1 (left or right), and a random waiting time 
Ti G [0, oo[. The total time i„, particle position a;„, and 
the spring extension Z„ follow the recurrence equations : 



tn — ^n— 1 ^" '^n 
•^n — -^n—l ^nd 
In — In-l—and + VTn 



(6) 



(4) 



A fundamental difference distinguishes constant force 
and spring friction experiments. In the first case, the 
waiting times are statistically independent, while in 
the second case, they are strongly correlated, because 
their probability distribution is conditioned by the spring 
extension l{t), which varies slowly. For example, a long 
waiting time extends the spring, and forces the next wait- 
ing times to be shorter than the average. These corre- 
lations prevent from performing a complete analytical 
study, and justify a numerical approach. By defining 
V = il3kdv)/2, X{t) = i(3kdxit))/2. Lit) = Xit) - Vt , 
A — (/3fcd^)/2, the cosine term becomes ch.{Vt + Lit)). A 
spring length L(t) ~ 1 means that the external force has 
reached its critical value Fc ~ 20. /d. L{t) is piece- wise 
linear, increasing during the pinning of the particle, and 
with discontinuities — A at each jump. A small value 
of A means that the particle has many different traps 
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accessible. At the opposite, for A > 1, the spring forces 
the particle to occupy only one given site or its nearest 
neighbours. This regime has not been considered here. 

Provided the temperature is not too low, ujq , uf^ , 
and the caracteristic time for the effective barrier are 
of the same order and constitute a unique, microscopic, 
time scale. We have restrict ourselves to this situation 
(/3 ~ 1). The only independent parameters of the model 
turn to be : those defining 7i, the rescaled jump length 
A and the "velocity" VcJ^^. 

The mean friction force is given by 



2T 

{F) = — {L{t)) 



(7) 



where brackets stand for the time average. 

(O) = lime_,oo ^ /o^ 0{s) ds. The exact formulas for the 

mean position L and its quadratic fluctuation are : 



{L) 



lim 

n— >oo 



lim 

n— i-oo 



j:tiL{t})'n + mVr! + V'ri/3 

En 



There is no need to average over a quenched disorder 
in this model, and simulations are notably simplified. We 
stress however that the time averages mentioned above 
may present important fluctuations from one realization 
of the process to another. This is especially the case when 
a broad distribution of trapping time is considered, and 
the displacement x{t) as a function of time may exhibit 
fluctuations as large as {x{t)) itself. Note that our nu- 
merical work corresponds to spring pulling experiments, 
thus trapping time distributions are regularized. We have 
checked carefully that the time averages were well defined 
and converged. 

At each jump, the direction cr„ = 1 is chosen with 
probability p(r„,t„_i) and Gn = —1 with probability 
1 — p. We made the choice 



p(T„,t„-l) 



exp(-/3/i+(r„)) 



exp(-/?/i+(T„)) + exp(-/3/i-(T„)) 
1 

"l + exp[-2L(t„_i)-2y -tJ' 



(8) 



Restricted to the case V — and H constant, the above 
choice fulfills the detailed balance equation, and the re- 
sulting Boltzmann equilibrium distribution for l{t) is a 
gaussian with variance T/k. 

One can demonstrate that the general case = 0, but 
Ti, increasing, still leads to a stationary distribution pro- 
file for l{t), provided the mean trapping time exists, i.e. 
Ti must be a bounded or slowly increasing function of 
time. 



IV. LARGE ALGEBRAIC WAITING TIME 
DISTRIBUTIONS. 

Definitions. It is well known that ageing effects ap- 
pear when the waiting time distribution decays like 
tt{t) ^ T--(i+") and a < 1. The following choice 
(ti = 1/uii is the reference time scale): 



1 



aTi 



, nit) 



1 

n 



1 



aTi 



-ia+l) 



is obtained, following (||), with 



Pn{t) = In ( 1 + 
Wtotir) 



and 



an -h T 

When a constant force F is applied, one has 



Wtot ( T ) 



ch(/3Fd/2). 



(9) 

(10) 
(11) 

(12) 



The effective exponent is no longer a, but oi!(F) — 
a ch(/3Fd/2) and the hopping length d remains constant. 
This picture is different from the one emerging from 
the Sinai model, where the force dependent exponent is 
a{F) = F/Fc {Fc is a given critical force), and the typi- 
cal hopping length depends on F like d{F) ^ 1/F^. This 
latter phenomenology could be used as another starting 
point instead of equations (^,§,||). 

The behaviour of a system with 71 growing as above, 
and a weak constant force F is just a particular case of 
the "continuous time random walks" described in [p^ . 
The mean displacement {x{t)) of the particle turns out 
to be in this case (disregarding prefactors): 



Fd 

{xiO) =«(()■<;. till — 



• F 



(13) 
(14) 



N(t) is the typical number of jumps during the time in- 
terval [0, i], and functions of /3Fd/2 have been linearized 
in the last equation. The sublinear dependence in time 
means that the velocity is asymptotically 0. 

Let us turn now to the system with a spring pulling 
the particle at constant velocity v. The constant force 
result (|lj) is helpless in this case. This is precisely a 
situation where driving with a spring is inequivalent to 
driving with a constant force. In this case, the broad 
waiting time distribution is cut when the spring depins 
the particle. A cut-off time tc is defined as : 



kl{tc) — kvtc — Fc — 



n{tc 



d/2 



ln(l + ic/aTi) kl3d 



(15) 



(16) 
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A mean field approach. One can try to understand 
the creep behaviour of the system by a mean field ar- 
gument, by neglecting all the fluctuations in the force 
and writing F ~ k{l(t)). Owing to the regularization 
provided by the velocity, the mean value of the waiting 
times exists and equals: 



{r{v)) 



Tl 



Tl 



(17) 



Result (13) holds, and one can write a self-consistent 
equation for v, v 1. 



2T t,{v) V n 



(18) 



which, apart from a logarithmic factor, leads to 
tc{v) ~ 1/v and : 



f3d^ \l3kdTi) 



(19) 



The linearized case (3Fd/2 < 1 implies a' = 
a ch(/3Fd/2) ~ a (at weak velocity, the force does not 
change much the value of a'). This mean field approach 
gives a power law creep F ~ i;". Unfortunately, it 
strongly disagrees with numerical simulations. 

The simulations were done, following the scheme of 
equation (^, section II. Details on the computational 
scheme are deferred to the Appendix : the main diffi- 
culty comes from the conversion of random numbers in 
waiting times, and a trick is necessary to avoid computer 
time wasting. 

Analysis of the numerical results. Figure (|^) shows 
an histogram for the distribution of waiting times with 
a = 0.5, at various velocities. The distributions are 
clearly cut at tc{v). The constant slope —{1 + a') in 
the log-log plot shows that one has true power laws until 
the cut-off. At intermediate velocities, a' departs from 
its velocity value, as expected. 

Figure (Q) shows the dependence of the v — F curve 
in the exponent a. The characteristics have been found 
to be monotonic and increasing. For both a = 0.9 and 
a = 1.3, the curves are presented with horizontal error 
bars (which are very small for a = 0.5 and a = 0.7). 
They give an estimate for the fluctuations between suc- 
cessive runs. The convergence is all the more difficult to 
get that the velocity is weak and the exponent a close 
to 1. All the characteristics with a < 1 have a thresh- 
old value decreasing while a becomes closer to 1. For 
a > 1, one should recover an ohmic regime at low veloci- 
ties. The characteristics for a = 1.3 are compatible with 
such an ohmic response, but fluctuations from one run 
to another remain important. This is probably a conse- 
quence of the infinite mean squared value of the waiting 
times when v —^ 

In log-log coordinates, the v — F curves are concave, 
and seems at first glance to end at a threshold value 



F{v 0) whereas in linear coordinates, at the opposite, 
they are convex. Among the curves known to change 
their convexity when the axes are logarithmically gradu- 
ated, one finds for example the celebrated stretched ex- 
ponentials V — vq exp[— (/c//)^] occurring in some glassy 
creep cases For a = 0.5 such a functional form 

fits well the intermediate regime 10"'' <v< IQ-i with 
7 ^ 2.8 ± 0.1, but does not describe anymore the slower 
regime v < 10~^. This latter case is better fitted by a 
threshold characteristics Fc + (F — F^,)^ , (3 = 2.1 (over 
only 1 decade). 

One can roughly say that the v — F curve, for a = 0.5, 
has a threshold value Fc with a force F which does not 
vary by more than 12 % over 2 decades. This is enough 
to discard the mean field description mentioned above 
(eq ^) . The other values of a = 0.7; 0.9 exhibits a sim- 
ilar threshold force, but as a becomes closer to 1, it is 
more and more difficult to get a good, well converged 
value for the force at weak velocity and the fit with the 
stretched exponential becomes poor. 

One can account for the existence of the quasi- 
threshold behaviour of the v — F curve {i.e. up to log- 
arithmic factors) with a simple argument. The distribu- 
tion of waiting times is a broad power law distribution 
cut for times greater than ~ d/v. The contribution of 
large waiting times t ^ d/v represents a finite ratio of 
the total sum t = X]r=i''"«- Such waiting times corre- 
spond to a finite value for the spring extension kvr, and 
occur at a constant rate. As they are not balanced by 
any negative contribution, the resulting mean spring ex- 
tension remains finite. Figure (^) confirms this scenario. 
The histogram for the variable h remains asymmetric, 
even for a vanishing velocity, leading to a constant fric- 
tion force. The "mean- field" approach fails because it 
relies on the existence of a typical value F for the force, 
vanishing with in contradiction with figure (^. 

It is interesting to notice that the mean-field procedure 
leads to a (wrong) power law creep at low velocity. Such 
power laws have been found by Horner in the mean-field 
theory of a particle in a short range correlated disor- 
der ijl^. Whereas it is a rather different problem, where 
ageing comes from slow relaxation in a high dimensional 
phase space, and not from any barrier mechanism (sec- 
tion U), we emphasize that both mean- field approaches 
lead to a power law for w — > 0. One can suspect the 
corrective terms to the above mean-field theory (which 
contain all thermal activation effects) to modify the al- 
gebraic creep into a threshold, like here. 

Anyway, in every physical case, it must exist a huge 
but finite equilibration time teq, at which the effective 
barrier TL stops increasing. The above threshold char- 
acteristics must match with a slow thermally activated 
ohmic regime for velocities v < d/teq. 
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V. VELOCITY INDUCED FLUCTUATIONS. 

A stepwise potential. The necessary and sufficient 
condition for the existence of an ohmic regime at low ve- 
locity is that the mean trapping time (r) remains finite. 
This is because the effective "viscosity" ?/ = lim^^o F/v 
scales like (r). At the opposite, (r) = oo implies a sub- 
linear or "creep" regime in the v — F characteristics. This 
latter can always be written: 



F — rj ■ V, 



(20) 



where 77 depends on or on w according to the friction 
experiment considered. When using a time-dependent 
potential Ti, one can generate a distribution of waiting 
times strongly dependent on with consequently, a non 
constant mean pinning time (r(u)), decreasing with v. 

Whenever TL remains bounded, one expects a crossover 
between the strong velocity regime involving 'H{t ~ 0) 
and the low velocity regime depending on the whole range 
T G [0, oo[. Moreover, a detailed study of the fluctuations 
{l'^{t)) of the particle's position turns out to be very in- 
structive and reveals some unexpected features. 

In order to make apparent the effects of a non-constant 
Ti function, I have introduced a decreasing step function. 
Such a sharp behaviour for Ti should make the crossover 
well contrasted. This is done with the escape rate : 



\ UJl forTe[0,ts] 
y U02 for r G [is, oo[ ' 



(21) 



where t^ is the arbitrary position of the step. The n(t) 
function (^ and its reciprocal 11" ^ can be computed ex- 
actly, providing an easy numerical generation of waiting 
times. n(t) equals : 



t <ts exp 



t > ts : exp 



V 
V 



(sh(L) - sh(L + Vt)) 

ijJ2 - OJl 



sh(L) 



V 



sh(X -f Vts) 



-ysh(L + l/t) 



(22) 



The enhancement of fluctuations. The mean square 
value is computed and compared with the equilib- 
rium one, given by the energy equipartition theorem 
^{^'^i't)) — where k denotes the spring stiffness. 

In figure (||) the ratio f = k(^P{t)') /T is plotted for 
the special choices uji — 1, uj2 — 0.01(= 1% oji), tg rang- 
ing from 0.5 to 7. The results are compared with the 
constant H cases A: w^^^ = oji and B: w^^^ = L02. The 
step function case S interpolates between these two static 
pinning potential cases. Thus, one expects the results for 
S to fall in between the two extreme cases A and B. 

The v—F characteristics behaves in this way (see figure 
(|)). Wher eas the constant w^^^ cases A and B exhibits 
fluctuations very close to equilibrium ones, regardless to 
the value w, the stepwise situation S shows a maximum 



for an intermediate regime centered around . Fluctu- 
ations are enhanced with a factor 3 or 4, before reaching 
their near equilibrium value at high velocities (figure |^). 

The value ~ 1 at high velocities suggests a decou- 
pling between translation motion and fluctuations, simi- 
lar to the one due to the galilean invariance, in the dis- 
order less case. This is reminiscent from a well known re- 
sult about a fast particle moving on a random potential. 
When considering a particle driven at high velocity, with 
a force greater than the critical one, the pinning poten- 
tial plays only a perturbative role. The random pinning 
force acts like a Langevin noise and leads to an effective 
"shaking temperature" scaling like the inverse velocity 
\/v (as far as a single particle is considered) . 

The effect of fluctuations is thus maximal at interme- 
diate velocities. In the context of a friction experiment, 
such anomalous fluctuations are reminiscent of "stick- 
slip" phenomena. On the other hand, in the physical 
case where the spring potential mimics the elastic inter- 
actions between particles, a strong fluctuation regime in- 
dicates that the motion occurs with strong deformations, 
with a possible breakdown of the elastic regime towards 
a plastic one. At the opposite, at weak or large veloci- 
ties, the motion of the particles may be coherent, within 
an homogeneous flow |17| (two neighbouring particles re- 
main close for a long time) . The onset of such a "plastic 
flow" should occur when the amount of fluctuations (x^) 
reaches a critical value like in the celebrated Lindemann 
criterion for melting. 

Figure(^) shows histograms for weak v = 10^^, inter- 
mediate V — Q.b and high velocity v = 2. The step occurs 
at is = 1- — 1 corresponds to a gaussian shape of the 
l{t) distribution, (p > 1 coincides with an asymmetric 
distribution for l{t), with a tail on the right, related to 
large values of I. The presence of such strong fluctua- 
tions requires a high enough velocity, able to "convert" 
a long pinning time into a large value of Z. At the other 
extreme, if v is too large, the friction force is higher than 
the critical force, and the pinning is no more efficient. 
The effect is maximal for 0.32 < vtg < 0.35. 



VI. CONCLUSION 

In this paper, a general stochastic process has been 
defined, which allows dynamical pinning effects, by in- 
troducing time-dependent barriers. The model has first 
been specialized to generate power-law distributed wait- 
ing times. The case of a particle driven by a spring 
has been investigated and compared to the previously 
known constant driving force case. The sublinear drift 
of the latter case is found to have a counterpart which 
in the present study appears as a finite threshold for 
the friction force, even for vanishing velocities. A naive 
mean-field approach fails in describing the velocity-force 
characteristics, showing the prominent role of large wait- 
ing times, even within a distribution cut-off by the drift 



6 



motion. The analogy with the D = (x "mean field" result 
suggests that the finite dimension corrections could drive 
the power-law creep into a threshold characteristics, ul- 
timately linear at very low velocities if the distribution 
of waiting times is cut by size effects. 

Then, we have investigated the consequences of a time 
dependent pinning potential concerning the fluctuations 
in the particle position. It has been achieved with a sharp 
step profile for the time dependent barrier. Fluctuations 
are strongly enhanced at an intermediate velocity regime. 
This phenomenon is important in various physical cases 
where the simple model presented here applies and de- 
serves more investigations. 

More generally, this model provides a general frame- 
work in order to test the equivalence between applied 
force and applied velocity in a friction experiment, with 
presence of "surface ageing" and at a scale where the 
temperature is a relevant parameter. 



has been tabulated, by a direct numerical calcu- 
lation, for a grid in the plane (L,liiV) with variable 
L between —3 < L < 5 ; AL — 0.1 and variable V, 
10~^ < y < 1 with 20 values following a geometrical 
sequence. At each point of the grid 11^ ^(r) is approx- 
imated with a Tchebitchev polynomials approximation, 
by keeping the 30 first coefficients. 

When a couple {L, In V) not on the grid is required, its 
Tchebitchev coefficients are linearly interpolated. The 
conversion r(r) is reduced to a fast polynomial evalua- 
tion. 

The two previous approximations match fine for 
Vt ~ 0.01, and were actually used in the simulation. 
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APPENDIX:RANDOM GENERATION OF 
WAITING TIMES. 

Given a random number r taken from a uniform dis- 
tribution over [0, 1], the corresponding waiting time r is 
n-i(r), with n-i(n(s)) = s and n defined by (|). As 
the integration cannot be -in the general case- performed 
analytically, a direct evaluation of r requires to find nu- 
merically the root of a function defined by a quadrature. 
Independently, convergence requires many millions (10^ 
to 10^) numbers at low velocities, excluding any direct 
computation scheme for 11^^. 

It's natural to look for asymptotic approximations of 
the integral. First, one notices that the integrand 
w^g^{T)ch.{Vt + L) has at least two very different time 
scales uj^^ and V^^. 

Let us detail the particular case ach{Vt + L)/{aTi+t). 
On the one hand, ii Vt ^ 1, the integral becomes 
ch(L) ads/{aTi + s), leading to : 



-(r) 



aTi 



exp 



ln(r) 
a ch{L) 



1 



(23) 



On the other hand, Vt > 1 describes situations where 
the force kl{t) approaches the critical depinning value. 
It's very unlikely to find high value L ^ 1 or nega- 
tive L <C — 1. In order to cover the range W ~ 1, 
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FIG. 1. Schematic illustration of the friction experiment. 
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FIG. 2. Barriers h+ and (Section 
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Histogram of waiting times 




log ( waiting time) 



FIG. 3. Distribution of waiting times regularized by the velocity, (a = 0.5). This curve H{\nt) has to be understood as 
follows : the number of waiting times lying between ti and t2 is N{ti,t2) = A/" ^ H{\nt)d\Tit. [Af is a normalization factor). 
H{lnt) ~ t-". 



V-F characteristics 
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FIG. 4. Dependence of the v-F curve in the exponent a. The closer to 1 is a, the more difficult is the convergence at low 
velocities. A = 0.005; V = 0.5f3kd ranges from 4.10"* to 0.1 
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Histogram for the spring length (l(t)) distribution. 



The area under the curve is 1 . 




-2.0 -1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 
Spring extension i(t). 

FIG. 5. The histogram of spring lengths L{t) (= (3kdl{t)/2) does not center around 0, even for a vanishing velocity (here 
V = 10~®). The inset shows the histogram with a logarithmic vertical axe. 
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FIG. 6. (j) = k(5 {l^{t)') ratio between fluctuations and equilibrium fluctuations. The stepwise H cases S lead to a mgiximum 
at t; = Vmax- The constant Ti cases A,B axe characterized by </!) ~ 1 (not on this picture). The inset shows the product Umax • ts 
vs ts, which has a constant value ~ 0.33. 
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FIG. 7. Comparison of the velocity force characteristics, between the constant w cases A and B, and the stepwise w cases S. 
The time ts of the step ranges from 0.5 to 5. 
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FIG. 8. The distribution of spring lengths L{t) is a gaussian at low and high velocities and is stretched at intermediate 
velocities. 
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